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Abstract 

We have presented an elegant high energy quantum problem, namely, the full Dirac oscillator under axial 
magnetic field with its full solution. We have found the energy spectrum which is rich and at the same 
time has a novel structure. The quantized energy levels show coupling of the oscillator frequency with the 
Larmor frequency in the 2D surface where the electrons under consideration follow a 2D oscillator. The axis 
in which magnetic field is pointed, the electrons follow a ID oscillator. There is also coupling between spin 
and orbital motion and also a coupling between a resultant effect of orbital and spin motion with Larmor 
precession. 
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1 Introduction 

Moshinsky and Szczepaniak [1] have introduced an interesting interaction in the Dirac equation such that the 
Dirac Hamiltonian is linear in both momentum and the coordinates. The resulting Dirac equation is known as 
Dirac oscillator because it turns into a harmonic oscillator with strong spin-orbit coupling in the non-relativistic 
limit. The oscillator is achieved by introducing the coupling in the Dirac equation as [1] 

ih-— = \cli. (y — + /3mc^l If", 

at 

where W is the four component bispinor corresponding to spin-half particle, like electron, c is the speed of light, 
and /3 are standard representation of Dirac matrices, w is the classical frequency of the oscillator. When 
w = 0 , we recover the ordinary Dirac equation[2]. 

The Dirac oscillator plays a vital role in the depiction of relativistic many body problems [3, 4 and 5] and 
super symmetric relativistic quantum mechanics [6, 7, and 8]. Dirac oscillator framework is also proposed in 
quantum chromodynamics, mainly in connection with quark confinement models in baryons and mesons [9]. 
Most thrilling aspects of this oscillator are its association to quantum optics [10]. It also maps into Anti-Jaynes- 
Cummiiigs model for describing the interaction of a atomic transition in two level systems [11, 12, and 13]. 
Moreover, a lot of attention and many physical applications in various branches of physics of this oscillator 
have been found, for example, in semiconductor physics [14], graphene [15] etc. Recently an experimental 
demonstration of this oscillator has been initiated [16]. 

In this context, we consider (3+1) dimensional Dirac oscillator in an external homogeneous magnetic field. 
We investigated mathematical formulation of the oscillator problem and solved it. The quantized energy spec¬ 
trum for up and down spin is derived. In the non-relativistic limit, this energy spectrum reduces to the spectrum 
of three dimensional harmonic oscillator with a spin-orbit contribution, when the magnetic field is absent. We 
present this problem and give a solution in section 2 and in section 3, we give a brief conclusion of the work. 


2 Dirac oscillator in magnetic field 


In presence of an external magnetic field, the (3+1) dimensional Dirac oscillator is given by 




^ + j3m(^ ip = Eip 


( 1 ) 
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where A is vector potential and —e(e > 0) is the charge of the oscillator which is considered to be an electron 
here. We consider the uniform magnetic field to be along the z- direction and vector potential to be in the 
symmetric gauge as ^ = (—"fy, -fx, O), where B is the magnetic field strength. Using the following standard 
form of lA and /3: 


and 




0 \ 
^ 0 ) 



and two-component form for ip as 



we obtain the following two simultaneous coupled equations: 


'ipb + imcojl^ 


and 


-I- Ipa — imcuj~^.l^lpa 

The system gives the following two equations 


[E - mc^) xpa 
[E + mc^) tpb- 


( 2 ) 

(3) 

(4) 

(5) 


[c^ + ce + ce (^•^) + 


imceoj 


^ -I- rn?(Auj‘^ ( = {E'^ — m^c^) -ipa 


and 


[c^ + ce + ce (^•^) + 


— imceuj 

One can easily show that 


( 0^.7^), ^ -I- ( o^.7^)^]V'6 = {E"^ — m^c^) ipb- 


( 6 ) 


(7) 


(^■^) = be, + Bha„ 




= A^ 


where electron spin ^ = | o^, 


and 




( 0^.1^), ^ = ir"^.1^ — .f) .f^] 

=r2. 


Using these, Eqs. (6) and (7) can be cast into 


-I- ceB{Lz + h<Tz) + e^A^ — — mceojr'^[azB — (o^.r) .f^] 


i’a 


and 


= — tAc^ + imc^fkv) Ipa 

Amc^uj - 


Ap^ + vAAAA + ceB{Lz + hoz) + e^A^ H —+ mceojA[azB — {~A.f) 


V’fc 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


( 12 ) 

(13) 


( 14 ) 
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= — Smc^huj) 'tfjb- (15) 

Since 1^ is along the z-direction, and in magnetic field spin of an electron processes aronnd 1^, the effective 
component of is dz. Hence, the second term within the square bracket in the 6 th term of the both equations 
should reduce as follows: 

(o^.f) = cTzB cos^ 9, 

where 9 is the azimuthal angle. 

Then we have 


22 | 2222 | TD / T 

c p +m c oj r + ceB[Lz 


fvJz) 


e^B^ 


sin^ 9 — 


Amc^u} - 




mceujazBr^ sin^ 6 


i’a 


and 


= + Smc^huj) i/'a 


99 9999 / B 9 9 ^ 4?nC UJ 

c p + mcLor + ceB{Lz + huz) H- —r sm 9 H- - — 

= [E'^ — — Smc^hu) ipb- 


tJ 


+ mceujCFzBr^ sin 6 


V’b 


(16) 


(17) 


The first two terms of the left hand side of Eqs.(16) and (17) constitute Hamiltonian of 3 dinensoinal harmonic 
oscillator. Third term reveals an interaction between 2 : component of total magnetic moment and the magnetic 
field and the fourth term reveals interaction between electric charge and the magnetic field. The fifth term 
appears due to spin-orbit interaction and sixth term due to an interaction between electron spin and magnetic 
field. 

We can find closed form solution of Eqs.(16) and (17) using either | Yj^’-{9^(f)xms > basis or | > 

basis, where the quantum number mi is associated with the eigenfunctions e*"*'*^ of Lz and SzXm, = fnshxmsj 
ms = ± 5 . In both of those bases, the eigen value of L is h^mimg. We prefer now to work with the 
I > basis. 

Since the system possesses cylindrical symmetry, in cylindrical coordinates p, (p and z , where +y^, 

Eqs. (16) and (17) can be written respectively as 


pdp \ dp 


2f,2 ) ^ d ( d 


and 


Amc^u! - 


1 d^ 9^1 2222 2222 ^ N 2 Amc^oj-p-^ 

r + rn P + m c to z + ceB(Lz + huz) H- -—p -r— L . S 

p^ dp dz^ j An 


— mceujazBp^Jipa = (E^ — m^c^ + ‘Hmc^hJ) pa 


ri d 
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\pdp 

ydp) 

1 _ 

9^2 

+ 


„2, 2 2 I ..^2 2, ,2.^2 


e'^B'^ 9 Amc^to 


(18) 


+ mceu}azBp^]pb = {E^ — mPc^ — imcPhw) pb- (19) 

As the equations involve partial derivatives of p and z and operators of z-component of angular momentum and 
spin, we take the solutions in the form : 

Pa = Eaip)Ga{z)E^‘'^Xm. = -l (20) 

and 


Pb = Fb{p)Gbiz)P^‘^Xm.=A- 


( 21 ) 


Then we get from Eq.(18) 

1 „2fc2l d (^dEa{p)\ ,_^2^2W 


c ^ -777 P 


Fa{p) pdp\ dp 


Let 


2„2 


B^e 


+ mceujB + mPcPuj"^ p^ — c^tP 


1 d^Ga{z) 
Ga{z) dz^ 


= E^ — m^m + imc^Huj — ceB{mi — l)h — 2mmujmih. 


E^ — m^c^ + Smc^huj — ceB{mi — l)h — 2mc^ujmih = A, 
where A is a constant. Hence, the Eq.(22) demands 


— c h' 


,2t2 1 9pGa{z) 2 2,2^2 


Ga{z) dz"^ 


+ m c UJ z = e, 


+ m^(?uj^z^ 

( 22 ) 

(23) 

(24) 
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where e is another constant. Then we have 


— c h' 


212 1 f (p) 


+ C^h^^Fa{p) + 
P 




p dp \ dp 

Multiplying Eq.(24) by 2 ;^, we have 

(PGa{z) 1 22 


+ mceujB + I p^Fa{p) + (e — X)Fa{p) = 0. (25) 


2m + 2 ““ " = W 

This is the simple harmonic oscillator Schrogringer equation and hence, 

e = 2 + -^ hjjjmc^, n = 0,1, 2... 

and Ga{z) is the wave function of simple harmonic oscillator. 

Eq.(25) can be written as 


Ga{z) 


(26) 


(27) 


<PFa{p) , 1 dFa{p) mf 

* ~p~zir ~ 7 ^ 


-( 


B'^e^ meijjB 


V Ac^h? IFc 




1 


P Fa{p) + ~ ^)Fa{p) — 0 


(28) 


Now, performing a change of variable 




^ 2 g 2 rneujB rrFuj'^ 




+ 


h?c 




-p, 


we obtain from Eq.(28) 


where 


d^FaiO , IdFaiO 


de 




D = 


-^Fa{0-eFa{0+DFa{0=^, 




A- 


^2fe2 / I mecuB i m?u)'^ 

^ " ''47^ + ~ 1F ^ + ~ 1 ^ 


According to [17], the corresponding solutions are 


(29) 

(30) 

(31) 


D = 2{N +l),N = 0,1,2... 


(32) 

and 

i"a(e)=Ae-HV™V,(0, 

where A is a normalization constant and /™(_|_j,(^) satisfy the equation 


(33) 

^d^fZl+kiO . , , dfZZZO . ,, 

^ + (Wi + 1 0 + kfmi+k 

o' 

II 

(34) 

. for N zs an 

2^ jo, 1,2. for N is an 

even integer 
odd intiger 

(35) 


and mi can take (iV/2) + 1 values for an even intiger N and A^ + 1 values for an odd intiger N. 

Equation (32) gives the energy spectrum of the (3+1) dimensional Dirac oscillator in presence of magnetic 
field: 


1 . 


= 2m<? (N + l)h{u; + ujl) + (n + ^)^ + ( ~ ^ fiw + {mi — 


for ms = -\ , 


where the Larmor frequency is 


UJL = 


eB 

2mc 


(36) 

(37) 


The Larmor frequency occurs in the energy spectrum due to interaction between electric charge and the magnetic 
field. The third term in the square bracket of Eq.(36) appears for spin-orbit coupling and the fourth term for 
an interaction between total magnetic moment and the magnetic field. 

Now, the other component of ip is 


V'fe = F,{p)G,{z)F^‘K^i = F,{0G,{z)F^^^Xm.=h = C{,^{0Gi,{z)F'^^’f> 


Xm. = i 


(38) 
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and the corresponding energy spectrum is 


E'^ — = 2m<? I (A^ + 1) (w + ujl) + (n + ^ fiw + {mi + l)hLUL 


for ms = \ . (39) 


For B = 0 


E"^ — = 2mc^ HN + n+ ^]huj+{mi — ^]haj 


= 2m(? 


3 

n' + -]huj+{mi-^]huj 


and 


£^2 _ ^2^4 ^ 2mc2 


for ms = 


N + n-\--]fujj+\mi-\--]fvjj 


(40) 


= 2m(? 


3 

n' + -\ fkij+[mi + -\ huj 


for ms = i , 


where 


n' = + n = 0,1, 2.... 


(41) 


(42) 


These (Eqs.(40) and (41)) are of the same nature as in spherically symmetric Dirac oscillator [18]; the difference 
occurred only because of the cylindrical symmetry we have assumed from the very beginning of the solution 
procedure. 

The non relativistic limit of the energy spectrum is obtained by setting E = m(? + K with the consideration 
K <1; mc^ and we get 


and 


For B = 0 


and 


K = {N + l)h{uj + ujl) + (n + 2 + {^i ~ ^ ~ = -^ 

K = {N + l)h{oj + ojl) + (n + 2)^ + ^ huj + {mi + l)hajL for ma = \. 

K = + ^ hw + - ^hw for ms = 

K = Sw + fiw for ms = \. 


(43) 

(44) 

(45) 

(46) 


Except spin-orbit contribution terms (mi — |) fiw in Eq.(45) and (mi -I- |) fiw in Eq.(46), these are exactly the 
energy spectrum of 3D non-relativistic harmonic oscillator. 

If we take the solutions in the form : 


i^a = Fa{p)Ga{z)E^‘'^Xm. = k 


and 


the corresponding energy spectrum of ipa and %pb are given by respectively as 


(47) 


(48) 


E^ — mfc^ = 2mc 


and 


E"^ — m^c^ = 2mc' 


{N + l)h{uj - ojl) + (n -I- ^)hw -[mi + ^]hio + {mi + l)fwjL 


for ms = k (49) 


{N + l)h{uj — lol) + (n -I- -)huj — [mi — ^]huj + {mi — l)hujL 


for ms = -i 


(50) 
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3 Conclusion 


We have solved the (3+1) Dirac oscillator with magnetic field and found the energy spectrum. The system is 
cylindrically symmetric: there are basically two oscillators, one is a 2D oscillator in the x — y plane and the other 
is a ID oscillator in the z- direction. The energy spectrum shows this splitting exactly. Moreover, there appears 
an oscillation which is basically a Larmor precession, coupled (i) with the 2D oscillator in x — y plane and (ii) 
coupled with a net effect of orbital and spin angular momentum. The 2nd coupling is shown in the last term 
in Eqs.(36) and (39) or (49) and (50). There is also spin-orbit coupling coupled with the Dirac oscillator in the 
third term in Eqs.(36) and (39) or (49) and (50). The non-relativistic limit of the energy spectrum is rightly that 
of 3D harmonic oscillators plus spin-orbit coupling. Everything appeared smoothly in our calculation because 
of the initial decomposition of the problem by a cylindrical symmetry. In the absence of magnetic field, the 
problem shows characteristics of a 3D Dirac oscillator as can be verified with what is available in literature. The 
energy spectrum is novel and rich. This oscillator might find application in quantum optics. The system can 
easily be manipulated using photon beams or electromagnetic signal and made to absorb and emit radiation 
with energy hu or huj^ or /i (w + lol) ov h{uj — ujl). This is clear from the expression of the energy levels. 

In conclusion, a very attractive quantum system has been elucidated which is a 3D Dirac oscillator immersed 
in axial magnetic field. 
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